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ABSTRACT 

The Fermi motion parameter is the most important parameter of ACCMM 
model, and the value Pp ~ 0.3 GeV has been used without clear theoretical or 
experimental evidence. So, we attempted to extract the possible value for pp 
theoretically in the relativistic quark model using quantum mechanical variational 
method. We obtained Pp ~ 0.5 GeV, which is somewhat larger than 0.3, and we 
also derived the eigenvalue of Eb — 5.5 GeV, which is is reasonable agreement 
with rriB = 5.28 GeV. 
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In the minimal standard model CP violation is possible through the CKM 
mixing matrix of three families, and it is important to know whether the element 
Vy5 is non-zero or not accurately. Its knowledge is also necessary to check whether 
the unitarity triangle is closed or not [1]. However its experimental value is very 
poorly known until now and its better experimental information is urgently re- 
quired. At present, the only experimental method to study V^b has been through 
the end-point leptonic spectrum of the B-meson semileptonic decays of CLEO [2] 
and ARGUS [3], and their data indicate that Vub is non-zero. Recntly it has also 
been suggested [4] that the measurement of hadronic invariant mass distribution 
in the inclusive B X^^y^lu decays can be useful in extracting \Vy,b\ with better 
theoretical understandings. In future asymmetric B factories with vertex detector, 
this will offer an alternative way to select h ^ u transitions that is much more 
efficient than selecting the upper end of the lepton energy spectrum. 

The simplest model for the semileptonic B-decay is the spectator model which 
considers the decaying 6-quark in the B-meson as a free particle. The spectator 
model is usually used with the inclusion of perturbative QCD radiative corrections. 
The decay width of the process B Xqlv is given by 



T,(B x,M . r.(6 ,M ^ |vwi^(^)/(5)li - . (1) 



where niq is the mass of the g-quark decayed from 6-quark. The decay width of 
the spectator model depends on m|, therefore small difference of m^, would change 
the decay width significantly. 

Altarelli et al. [5] proposed their ACCMM model for the inclusive B-meson 
semileptonic decays. This model incorporates the bound state effect by treating 
the 6-quark as a vitual state particle, thus giving momentum dependence to the 
6-quark mass. The virtual state 6-quark mass W is given by 



in the B-meson rest frame, where rusp is the spectator quark mass, ttib the B-meson 
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mass, and p is the momentum of the 6-quark. 

For the momentum distribution of the virtual 6-quark, AltareUi et al. consid- 
ered the Fermi motion inside the B-meson as the Gaussian distribution 



with a free parameter of Gaussian width. And the decay width is given by 
integrating the width F^, in (1) with the weight 0(p). Then the leptonic spectrum 
of the B-meson semileptonic decay is given by 



where Pmax is the maximum kinematically allowed value of p = |p|. The ACCMM 
model, therefore, introduces a new parameter pp for the momentum measure of the 
Gaussian distribution inside B-meson instead of the 6-quark mass of the spectator 
model. In this way the ACCMM model incorporates the bound state effects and 
reduces the strong dependence on 6-quark mass in the decay width of the spectator 
model. 

The Fermi motion parameter Pj^ is the most essential parameter of the AC- 
CMM model as explained above. However, the experimental determinations of its 
value from the leptonic spectrum have been very ambiguous until now because the 
various effects from the input parameters and the perturbative QCD corrections 
are intermixed in the spectrum. The value Pp ~ 0.3 has been widely used in ex- 
perimental analyses without theoretical or experimental clean justification, even 
though there has been recently an assertion that the BSUV model [6] is approx- 
imately equal to ACCMM model at Pp ~ 0.3. Therefore, it is strongly required 
to determine the value of Pp more firmly when we think of the importance of its 
role in experimental analyses. The better determination of Pp is also interesting 
theoretically since it has the physical correspondence related to the Fermi motion 
inside B-meson. In this context we are going to theoretically determine the value of 
Pp in the relativistic quark model using quantum mechanical variational method. 
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(3) 




(4) 
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We consider the Gaussian probability distribution function 0(p) in (3) as the 
absolute square of the momentum space wave function x(p) of the bound state 
B-meson, i.e., 



?.(p) =4ir|x(p)r , (5) 

The Fourier transform of x(p) gives the coordinate space wave function which 
is also Gaussian, 

^(r) = (^)3/2e--V,/2 . (7) 

Then we can approach the determination of in the framework of quantum 
mechanics. We apply the variational method with the Hamiltonian operator 

= ^p2 + + ^p2 + m2 + V{r) (8) 



and the trial wave function 



^(r) = (^)3/2e-AV2 , (9) 



where is the variational parameter. The ground state is given by minimizing the 
expectation value of i?, 

^E) = {i,\H\i,) = E{^) , (10) 
-^E{ii) = at = , (11) 

and then — Pp and E = E{p,) approximates uib- The value of /i or pp. corre- 
sponds to the measure of the radius of the two body bound state as can be seen 



4 



from 

(r) = Ai , and (r^)^ = ■ (12) 

In (8), for simplicity, we take the Cornell potential which is composed of the 
coulomb and linear potentials, 

V{r) = -y + Kr . (13) 

For the values of the parameters ac (= |tts), K, and the 6-quark mass m^, we use 
the values given by Hagiwara et al. [7], 

ac = 0.47 (as = 0.35), K = 0.19 GeV'^, = 4.75 GeV^, (14) 

which have been determined by the best fit of the (cc) and {bb) bound states. We 
will also consider the following value of etc for comparison in our analysis 

ac = 0.32 {as = 0.24) , (15) 

which corresponds to as{Q^ = "^s)- 

Before applying our variational method with the Gaussian trial wave function 
to the B-meson system, let us check the method by considering the T{bb) system. 
The Hamiltonian of the T{bb) system can be approximated by the non-relativistic 
Hamiltonian 

H c^2mb+ — + V(r) ■ (16) 

rrib 

With the parameters in (14) or (15), our variational method with the Gaussian 
trial wave function gives Pp — U — 1.1 GeV and E = -E(/i) = 9.49 GeV . Here 
Pp = 1.1 GeV corresponds to the radius R{T) = 0.2 /m, and ^(T) = 9.49 GeV is 
within 0.3 % error compared with the experimental value E'exp = i^T — 9.46 GeV . 
Therefore, the variational method with the non-relativistic Hamiltonian gives the 
fairly accurate results for the T ground state. 
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However, since the u- or d- quark in the B-meson is very hght, the non- 
relativistic description can not be apphed to the B-meson system. For example, 
when we apply the variational method with the non-relativistic Hamiltonian to the 
B-meson, we get the following results 

= 0.29 GeV, E = 5.92 GeV for ag = 0.35, (17) 
= 0.29 GeV, E = 5.97 GeV for as = 0.24. (18) 

The above masses E are much larger compared to the experimental value uib — 
5.28 GeV, and moreover the expectation values of the higher terms in the non- 
relativistic pcrturbative expansion are bigger than those of the lower terms. There- 
fore, wc can not apply the variational method with the non-relativistic Hamiltonian 
to the B-meson system. 

We use the following Hamiltonian for the B-meson system in our analysis by 
treating u- or d-quark relativistically, 

H-M+^ + ^p^ + m^ + V{r) , (19) 

where M — mi, and m — m sp- 
in our variational method the trial wave function is Gaussian both in the 
coordinate space and in the momentum space, so the expectation value of H can 
be calculated in either space, 

(i/) = (V^(r)|i/|V^(r)) = (x(p)|i/|x(p)) ■ (20) 

Also, the Gaussian function is a smooth function and its derivative of any order 
is square integrable, thus any power of the Laplacian operator is a hermitian 
operator at least under Gaussian functions. Therefore, analyzing the Hamiltonian 
(19) with the variational method can be considered as reasonable even though 
solving the eigenvalue equation of the differential operator (19) may be confronted 
with the mathematical difficulties because of the square root operator in (19). 
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With the Gaussian trial wave function (6) or (9), the expectation value of H 
can easily be calculated besides the square root operator, 

(p') = mr)\p'mr)) = (x(p)Ip'Ix(p)) = , (21) 
{V{r)) = (V'(r)l - ^ + XrlV'(r)) = -^{-ae^l + K/y^) . (22) 



Now let us consider the expectation value of the square root operator in the mo- 
mentum space 

(V^p^Tm^) = (x(p)|Vp^Tm2|x(p)) 

oo 

^ (^)Y e-P'/^Vp2 + ^2 d^'p 



oo 

= e-^y^s + (m//i)2 x^dx . (23) 



The integral (23) can be given as a series expansion by the following procedure. 
First, define 

oo oo 

I{s) = J + s x^e-'^'dx ^s^ j Vt^ + 1 t^e-'^'dt , (24) 



oo oo 

(s) = J y/x'^ + s e-'^'dx J V^Tl e''^^ dt . (25) 



h{s) = 





Next, from (24) and (25), we find the following differential relations 

Combining (26) and (27), we get a second order differential equation for I{s), 

sl"{s)-{l + s)l'{s) + h{s) = . (28) 
The series solution to eq. (28) is given as 



I{s) = ci/i(s) + 02/2(5) , 
I,is)^s'F{s;l3)^s'{l + ls + ^s' + 



192 



1024 



•}' 



(29) 



-— s^lnsf 1 + + + • 



1 

16" 



32 



1 5 2 
-s H s 

2 32 



192 



1536 



where F(s; |,3) is the confluent hypergeometric function which is convergent for 
any flnite s, and the integral constants c\ ~ —0.095, C2 = —1. See Appendix for 
the derivation of these numerical values for Cj. 

Finally, collecting (21), (22) and (23), the expectation value of B. is written as 



+ 



2/i 



1 + 



\{^ll^f + (i - 2ci)(W/^)' + \{mli,f\^{mli,) 



(30) 



up to {m/ jj)^. 

With the input value of m = rusp = 0.15 GeV, we minimize {H) of (30), and 
we obtain 



fx^OMGeV, ruB^E = 5.54: GeV for = 0.35 , 



fl = 0.49 GeV, 



E = 5.63 GeV 



for a, = 0.24 . 



(31) 



For comparison, we calculated (H) for the case of m = in which the integral of 
the square root operator is exact, and we get 



p = 0.53 GeV, E = 5.52 GeV for as = 0.35 , 
jl = 0A8GeV, E = 5.60 GeV for = 0.24 . 



(32) 



The calculated values of the S-meson mass, E, are much larger than the measured 
value of 5.28. The large values for the mass are originated partly because the 



Hamiltonian (30) does not take care of the correct spin dependences for B and B*. 
The difference between the pseudoscalar meson and the vector meson is given arise 
to by the chromomagnetic hyperfine sphtting, which is 



3Mm 



^ V72/ "^Cx 

Si • S2 V ( ) . 



(33) 



Then the expectation value of Vg is given as 



aciJr 



v/tt Mm 
2 acjJ^ 



30? Mm 



for B , 
for B* , 



and we treat (Vs) only as a perturbation. And we get for B meson 



(34 - 1) 
(34 - 2) 



= 0.54 GeV, Eb = 5.42 GeV for as = 0.35 , 
= 0.49 GeV, Eb = 5.56 GeV for as = 0.24 . 



(35) 



The perturbative result for is 



= 0.54 GeV, Eb' = 5.58 GeV for as = 0.35 , 
p^ = 0.49 GeV, Eb' = 5.65 GeV for as = 0.24 . 



(36) 



The calculated values of the 5-mcson mass, 5.42GcV («s = 0.35) and 5.56GeV 
(ttg = 0.24) are in reasonable agreement compared to the experimental value of 
m_B = 5.28GeV; the relative errors are 2.7% and 5.3%, respectively. But for Fermi 
motion parameter p^, the calculated values, 0.54GeV {ag — 0.35) and 0.49GeV 
{as — 0.24), are somewhat larger than the value 0.3, widely used in the exper- 
imental analyses of energy spectrum of semileptonic B meson decay. The value 
Pp — 0.3 corresponds to the S-meson radius Rb ~ 0.66 fm, and it seems too large 
to us. On the other hand, the value = 0.5 corresponds to Rb ~ 0.39 fm, which 
looks in reasonable range. 
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The Fermi motion parameter pp is the most important parameter of ACCMM 
model, and the value Pp ~ 0.3 GeV has been widely used in experimental anal- 
yses without clear theoretical or experimental evidence. Therefore, it is strongly 
required to determine the value of Pp more firmly when we think of the importance 
of its role in experimental analyses. We attempted to extract the possible value 
for pp theoretically in the relativistic quark model using with quantum mechanical 
variational method. The better determination of Pp is also interesting theoreti- 
cally since it has the physical correspondence related to the Fermi motion inside 
B-meson. We obtained Pp ~ 0.5 GeV, which is somewhat larger than 0.3, and we 
also derived the ground state eigenvalue of Eb — 5.5 GeV, which is is reasonable 
agreement with tub — 5.28 GeV. 
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The integration constants ci and C2 in (29) are given by the following relations. 
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Appendix 



oo 




{A -I) 



I"{S 0) 





2ci + C2(--^ Ins 



32 
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oo 




{A -2) 
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Then, from (A-1), we get 



C2 = -1 . 



The integral in (A-2) can be expanded as 



(A -3) 



J[s = a^)= x\x^ + a^r"\-^'' dx 



oo 

j x'^[{x + af - 2ax\ 



-3/2 



e ^ dx 



x'^(x + a) ^ 



1 - 



2 ax 



1 -3/2 



{x + ay 



e ^ dx 



E 

n=0 



(2n + l)!a 



oo 



2"(n!)2 J (x + a)2«+3 



e ^ dx 



{A -A) 



Next the integral in (A-4) is obtained by 



oo 

I 



X 



n+2 



(x + a) 



2n+3 



e ^ = 



1 



■e ^ dx . 



(2n + 2)! \da J J x + a 



{A -5) 



Again the integral in (A-5) is related to another integral, for a small value of a. 



n k=0 ^ ^ n 



The integral in (A-6) can be expanded in a similar way as to obtain the series (29) 
by making use of differential eqations. For a small value of a, 



00 2 

/ ^dx = --e-«' 
J x + a 2 



(21na + 7 + a2 + i|- + ---) + 0Fe "'(a + + + • • •) , 



3" 5 2! 

{A -7) 

where 7 ~ 0.5772 is the Euler's constant. In this way the constant ci is given by 
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an infinite series, 



oo 



3 1 

ci = -TTT + 5: - ttV^ ~ -0.0975 . (^-8) 
64 16 8 ^n2^ ^ ^ 

n=l 
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